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Abstract 

With a densely defined symmetric semi-bounded operator of nonzero 
defect indexes Lq in a separable Hilbert space H we associate a topo- 
logical space Ol (wave spectrum) constructed from the reachable sets 
of a dynamical system governed by the equation uu + (Lq)*u = 0. 
Wave spectra of unitary equivalent operators are homeomorphic. 

In inverse problems, one needs to recover a Riemannian manifold 0, 
via dynamical or spectral boundary data. We show that for a generic 
class of manifolds, £1 is isometric to the wave spectrum $1l of the 
minimal Laplacian Lq = — ^\c°°(Ci\dO,) acting in % = L2(£l), whereas 
Lq is determined by the inverse data up to unitary equivalence. Hence, 
the manifold can be recovered (up to isometry) by the scheme 'data 

T „ isom „, 

=> l =>- n Lo = n\ 

The wave spectrum is relevant to a wide class of dynamical sys- 
tems, which describe the finite speed wave propagation processes. The 
paper elucidates the operator background of the boundary control 
method (Belishev'1986), which is an approach to inverse problems 
based on their relations to control theory. 



Introduction 



0.1 Motivation 

The paper introduces the notion of a wave spectrum of a symmetric semi- 
bounded operator in a Hilbert space. The impact comes from inverse prob- 
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lems of mathematical physics; the following is one of the motivating ques- 
tions. 

Let Q be a smooth compact Riemannian manifold with the boundary T, 
—A the (scalar) Laplace operator, Lq = — A\c°°(n\r) the minimal Laplacian 
in % = /^(fi). Assume that we are given with a unitary copy Lq = ULqII* 
in a space H = WH (but Q,7-L and U are unknown!). To what extent does 
Lq determine the manifold Q? 

So, we have no points, boundaries, tensors, etc, whereas the only thing 
given is an operator Lq in a Hilbert space %. Provided the operator is 
unitarily equivalent to Lq, is it possible to 'extract' Q from Lq! Such a 
question is an invariant version of various setups of dynamical and spectral 
inverse problems on manifolds [2J, [3]. 

0.2 Content 

Substantially, the answer is affirmative: for a generic class of manifolds, any 
unitary copy of the minimal Laplacian determines Q up to isometry (Theorem 
1). A wave spectrum is a construction that realizes the determination Lq =>- Q 
and thus solves inverse problems. In more detail, 

• With a closed densely defined symmetric semi-bounded operator Lq of 
nonzero defect indexes in a separable Hilbert space % we associate a topo- 
logical space Ql (its wave spectrum). The space consists of the atoms of a 
lattice with inflation determined by Lq. The lattice is composed of reachable 
sets of an abstract dynamical system with boundary control governed by the 
evolutionary equation u tt + LqU = 0. The wave spectrum is endowed with a 
relevant topology. 

Since the definition of Ql q is of invariant character, the spectra Ql and 
of unitarily equivalent operators Lq and Lo turn out to be homeomorphic. 
So, a wave spectrum is a (hopefully, new) unitary invariant of a symmetric 
semi-bounded operator. 

• A wide generic class of the so-called simple manifolds is introduced^- The 
central Theorem 1 establishes that for a simple Q, the wave spectrum of 
its minimal Laplacian Lq is metrizable and isometric to Q. Hence, any uni- 
tary copy Lq of Lq determines the simple Q up to isometry by the scheme 



Roughly speaking, a simplicity means that the symmetry group of Q is trivial. 
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L =>- f2j = Q Lo = fi. In applications, it is the procedure, which recov- 
ers manifolds by the boundary control method [2],jl]: concrete inverse data 
determine the relevant L , what enables one to realize the scheme. 

• We discuss one more option: elements of the space % can be realized as 
'functions' on VLl IE Hopefully, this observation can be driven at a functional 
model of a class of L s and/or Spaces of Boundary Values. Presumably, this 
model will be local, i.e., satisfying supp (L™ d )*y C suppy. 

0.3 Comments 

• The concept of wave spectrum summarizes rich 'experimental material' 
accumulated in inverse problems of mathematical physics in the framework 
of the BC-method, and elucidates its operator background. For the first 
time, Ql has appeared in [1] in connection with the M.Kac problem; its 
later version (called a wave model) is presented in [I] (sec 2.3.4). Owing to 
its invariant nature, Ql promises to be useful for further applications to 
unsolved inverse problems of elasticity theory, electrodynamics, graphs, etc. 

Our paper is of pronounced interdisciplinary character. 'Wave' terminol- 
ogy, which we use, is motivated by close relations to applications. 

• The path from Lq to Ql passes through an intermediate object, which 
is a sublattice of the lattice £("H) of subspaces of the space %. Section 1 
is an excursus to the lattice theory, in course of which we introduce lattices 
with inflation. The wave spectrum appears as a set of atoms of the relevant 
lattice with inflation determined by L . 

• We give attention to connections of our approach with C*-algebras. As 
is shown, if Q is a compact manifold then Ql is identical to the Gelfand 
spectrum of the algebra of continuous functions C(Q). By the recent trend 
in the BC-method, to recover unknown manifolds via boundary inverse data 
is to find spectra of relevant algebras determined by the data [5]. We hope 
for utility and further promotion of this trend. 

Acknowledgements I would like to thank A.B.Alexandrov for kind consul- 
tations, and V.I.Vasyunin and N.Wada for useful discussions on the subject 
of the paper. 



2 In the BC-method, such an option is interpreted as visualization of waves [4]. 
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1 Lattices with inflation 

Reducing the volume of the paper, we do not prove Propositions. The proofs 
are quite elementary and typical technique is demonstrated in Appendix. 

1.1 Basic objects 

1. Lattice. Let £ be a lattice, i.e. a partially ordered set (poset) with the 
order ^ and operations a A b = inf {a, b}, a V b = sup{a, b}. Also, we assume 
that £ is endowed with the least element satisfying < a for a ^ [TJ. 

The order topology on £ is introduced through the order convergence: 
Xj — > x if there are the nets {fljjjejt an d {bj}je.ii (J is a directed set) such 
that dj ^ Xj ^ bj and sup{cij} = x — inf{6j} holds (we write dj t % and 
bj \, x). For an A C £, the inclusion x E A occurg^l if and only if there are 
dj ,bj £ A such that dj j" x and /or 6j 4 x [7J . 

Example 1. The lattice £ = 2 n of subsets of a set f2 with the order 
^=^, operations A = D, V = U, and = 0. 

Example 2. The (sub)lattice O C 2 n of open sets of a topological space 
Q. The convergence Uj "| w means w = sup{wj} = UjWj. The convergence 
ujj i u means u = inf {cjj} = int Dj Uj, where hit A is the set of interior points 

of a c a 

2. Inflation. For a lattice £, the set J-£ := ^([0, oo); £) of £-valued func- 
tions is also a topologized lattice with respect to the point-wise order, oper- 
ations, and topology. 

Definition 1. A map /:£—)■ is said to be an inflation if for all 
a, b G £ and s, t G [0, oo) one has 

(i) (ia)(0) = a and J0 £ = Ojr 

(ii) a < 6 and s < i imply (Ia)(s) < (76) (t). 
Inflation is injective: I~ 1 f = /(0) on /£. 

Example 3. is a metric space with the distance d. For a subset A G Q, 
by A t := {x G Q} \ d(x, A) < t} (t > 0) we denote its metric neighborhood, 
ant put A := A, 0* = 0. The map M : 2 n -»■ J" 2 n, (MA)(t) := A*, i ^ is 
a metric inflation. The image M2 Q is a semilattice: Ma V M6 = M(a V b) G 
M2 Q . The image of open sets is a (sub) semilattice MO C C J-^n. 

3 Everywhere ( ) denotes a topological closure. In some places, to avoid the confusion, 
we specify the space. 
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3. Atoms. Basic construction. Let V be a poset with the least element 
0. An a E V is called an atom if < a ^ a implies a = a [7]. By AtV we 
denote the set of atoms. 

Example 4. Each atom of 2° is a single point set: At 2 n = {{x} \ x E fl}. 
Example 5. If the open sets of a topological space Q are infinitely di- 
visibkS then At O = 0. 

Inflation preserves atoms: I At £ C At /£. 

For any lattice with inflation, the set Qjc '■= At /£ C J-V (the closure 
in topology of is well defined 0. This set is a key object of the paper. 
Namely, the following effect will be exploited: there are lattices and inflations 
such that At £ = but At 12, ^ 0. Inflation can create atoms! 

There is a natural topology on Qjc C J- £ . For atoms a, (3 E flic, we say 
that a influences on (3 at the moment t if a(t) A (3(e) ^ 0^ for any positive 
s. Define t aP := inf{t ^ | a(t) A (3(e) ^ £ \/e > 0}. If a(t) A /3(e) = £ for 
all positive i and e, we put £ aj g = oo. 

A function Tj£ : Vtjc X Vt IC —> [0, oo) U {oo}, Tiz(a, (3) := max{t a p, t^} is 
called an interaction time. 

Define the 'balls' B r [a] := {/3 E Vt IC \ T I& (a,0) < r} (r > 0), B°[a] := a. 

Definition 2. By (fij£, rzx) we denote the topological space that is the 
set Qjc endowed with the minimal topology, which contains all balls. 

Surely, at this level of generality, to expect for rich properties of this space 
is hardly reasonable. However, in 'good' cases the function tjz turns out to 
be a metric. 

Proposition 1 Let (O, d) be a complete metric space, 2 = 2 Q , I = M (see 
Example 3). The correspondence fi9if> M{x} E J"2 n a bisection between 
the sets tt and Vt M2 n = At M2 n = At M2 n = MAt 2 n = {M{x} | x E f2}. 
The equality r M2 n(M{x} ) M{y}) = d(x,y) holds. Function r M2 n is a metric 
on atoms, whereas (n M2 f!,r M2 fi) is a metric space isometric to (O, d) . The 
isometry is realized by the bisection M{x} -H- x. 

There are another topologies on atoms, which are also inspired by the 
metric topology. The first one is introduced via closure operation: for a set 

4 i.e., for any there is ^ B e O such that B C A and A\B ^ 0. 

5 But the case fij£ = is not excluded. 
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W C fi/£, we put 



It is easy to check that the map W H- W satisfies all Kuratovsky's axioms 
and, hence, determines a unique topology pjc in Qjc- Note a certain resem- 
blance (duality) of such a topology to Jacobson's topology on the set X of 
primitive ideals of a C*-algebra A. Namely, for a W C X, one defines its 
closure by 

W:= ji GX 

(see, e.g., [15]). 

One more topology is the following. We define the 'balls' by 




3t = t Q (a, 0, r) > s.t. ^ /3(t ) I a(r) \ (r > 0). 



As one can verify, the system {B r [a}} a£ n IC ,r>o is a base of topology. Hence, 
it determines a unique topology that we denote by am- 

If C = 2 K " and / is the (Euclidean) metric inflation, the topologies 
tic-, Pic, and ajc coincide with the standard Euclidean metric topology in 



4. Isomorphic lattices. Let £ and £' be two lattices with inflations 
/ and /' respectively. We say them to be isomorphic through a bijection 
i : £ — > £' if i preserves the order, lattice operations, and i(IA) = I'i(A) 
holds for all A e £. 

The bijection % is extended to bijection on functions % : — > by 
(if)(t) := i(f(t)), t ^ 0. The following fact is quite obvious. 

Proposition 2 lattices with inflation £ and £' are isomorphic then the 

spaces (^lic, T is.) ond (fV/y, Tr'£') flre homeomorphic. The homeomorphism 
is realized by the bijection i on atoms. 



1.2 Lattices in metric space 

5. Lattice (9. Return to Example 3 and assume in addition that 
Al. Q is a complete metric space 
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A2. All the balls {x} 1 are compact and {x} l \{x} s ^ as s < t. 

By A2, open sets are infinitely divisible. Therefore, we have At O = 0. 

Fix am G H and define the functions x*, x* G To'- x*(t) := {x}' as t > 0, 
x*(0) := 0o, and x*(t) := int{x}* as t ^ 0. Evidently, we have x* ^ x* in 
To- The upper function satisfies x* = lim e ^ M ({x } £ ) G M O, x*(0) = 0©. 
The 'clearance' between the functions is small: x*(t) = x*(t), i ^ 0. 

Since x* g MO, the segment [x*,x*] := {/ G 7*o | x* < / ^ x*} inter- 
sects with MO. The poset [x*,x*] H MO is a closed subset in J-© bounded 
from below. Hence, it contains minimal elements, which can be easily recog- 
nized as the atoms of MO. So, Vt M o ■= AtMO ^ 0. 

Example 6. For f2 C R™ one has x* = x*. Therefore, each segment 
[x*, x*] contains one (and only one) atom {x}*, t ^ 0. We don't know whether 
the same is correct for a Riemannian manifold Q. 

For an atom a g AtMO, define a kernel a := n t>0 C f2. 

Proposition 3 For each a, its kernel a consists of a single point x a g fl. 
Each atom a belongs to the segment [(x a )*,(x a )*]. If a — $ then a(t) = 
J(tj, t ^ holds. 

These facts follow from a general lemma stated and proved in Appendix. 

With each x g Q one associates the class of atoms (a) x := [x*,x*] fl 
AtMO. For a, /3 g (a) x one has a(t) = f3(t) (= {x}' ), t ^ 0. Hence, a and 
/3 interact at any t > 0. As a result, we have Tuoip--, P) — 0. 

The relation {a ~ /?} {tmo(o;,/3) = 0} is an equivalence on f^MO- 
The factor-set fi' MO := J1mo/~ is bijective to f2 through the map (a) >->■ a. 
The function r' MO ({a), (/?)) := tmo(<x,P) is a metric on fi^ . The equal- 
ity T' MO ({a) x , (cx)y) = d(x,y) is valid for all x,y g f2 and we conclude the 
following. 

Proposition 4 TTie metric space (0' MO ,r^ o ) is isometric to (f2,d). JTie 
isometry is realized by the bijection (a) x -H- x. 



6. Lattice O reg . For a set A C denote by OA := AnQ\A its boundary. 
Note that d(A C\B) C dAUdB and d(A U B) CdAU dB. It is convenient 
to put dtt = 90 = 0. 

Recall that we deal with complete and locally compact metric spaces (see 
Al,2). In addition, assume that Q is endowed with a Borel measure \i such 
that 
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A3. For any icO and t > 0, the relation ^(dA 1 ) = holds. 

Example 7. Q is a smooth Riemannian manifold with the canonical mea- 
sure (volume). In particular, Q C M n with the Lebesgue measure [TT] . 

Definition 3. An open set A C C is called regular^ if /i(cM) = 0. The 
system of regular sets is denoted by C reg . 

As is evident, O reg is a sublattice in O. It is a base of (9: each open set 
is a sum of regular sets (balls). By A3, O rcg is invariant w.r.t. the metric 
inflation: (MC reg )(t) C O reg , t ^ 0. In other words, we have MO reg C J 7 ^. 



Fix an x G fi. Note that x*,x* G J*o«g and x* G MO reg . Using the 
arguments quite analogous to ones, which have led to Proposition HI and 
factorizing the set of atoms w.r.t. the same equivalence ~, one can arrive at 
the following result. 

Proposition5 The metric space (Q' MO r Cg ,T' MOTCg ) is isometric to (Q,d). 
The isometry is realized by the bijection (a) x <-> x. 

The operation A h- Y A* := int(f2\74) is well defined on O reg and called a 
pseudo- complement [7J. The relations A fl A* = and A C (A*)* are valid. 

7. Lattice Introduce an equivalence on O reg : we put A ~ B if A = 5. 
Define := rcg / ~ . By [A] we denote the equivalence class of A. 

Endow Dl with the order and operations: 
[A] ^ [B] if A C B 

[A] A [B] := [ARB], [A] V [5] := [A U B] 
[A]^ := [A*] (= [int(0\A)]). 

The least and greatest elements are := [0] and 1 := [fi]. 

One can easily check the well-posedness of these definitions and prove the 
following relations: 
[A] A = 0, [A]V[A] ± = 1 

([A] A [B]) 1 - = [A] 1 V [B]\ ([A] V [fi]) X = [A] 1 A [J3J-L. 
Hence 9i is a lattice with the complement [7J. 

For / G J- orcg , define [/] G 7^, [/](*) := [/(*)], t > 0. 

Introduce the metric inflation on d\ by M : £H — > J-~<y\, (M[A\) (t) := 
[{MA)(t)] = [A% t^0. ' 

The relation At Mil = {[a] | a G AtM0 re s} holds. The map A H> [A] 
identifies the atoms belonging to the same class: if a, (3 G then a;(i) = 



6 Our definition is similar to (but differs from) the definition of regularity in [7], p 2f6. 
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(3(t), t ^ that implies [a] = [/?]. By this, the set SImvk = AtM9t is bijective 
to fl, whereas the 'interaction time' tmsk turns out to be a metric. 

Proposition 6 The metric space (flMm, T Mm) is isometric to (fl,d). The 
isometry is realized by the bisection [a] -H- x a . 

8. Lattice Dl n . Introduce a Hilbert space H := L 2)M (fi). 

For a measurable set A C fi, define the subspace HA := {xaU \ y €.H}, 
where xa is the indicator of A. Such subspaces are called geometric. If 
A G O reg then n(A\A) = fi(dA) = that leads to HA = HA. 

Definition 4. If A G (9 reg , we say the subspace HA to be regular. The 
system of regular subspaces is denoted by Dl H . 

Let £(H) be the lattice of subspaces of the space H (see item 10 below). 
The system Dt^ C £(H) is a sublattice. 

Introduce a map z : 91 — > 9^, iL4] := HA. As is easy to check, it 
preserves the operations and complement^. 

Extend i to functions: for an / <G we put if G T^h C J-2,{H)i {if)(t) '■= 
i(f(t)), t ^ 0. Also, define a metric inflation on 9^ w by iM : 9l n — >■ J-z(h)i 
(iMHA) (t) := H ((MA) (t)) = i[A*] = HA\ t > 0. 

Thus, i is an isomorphism of lattices with inflation. Propositions [3] and [6] 
lead to the following result. 

Proposition 7 The metric space (QiM<x H i T iNm n ) is isometric to (fl, d). The 
isometry is realized by the bijection i[a] 4-> x a . 

The meaning of the passage C reg — > 9l n is that it 'codes' regular sets 
in Hilbert terms. Later in inverse problems, we will determine the Hilbert 
lattices from inverse data, and then 'decode' them, i.e., extract information 
about geometry of fl. 

9. Dense sublattice. We say a system of subsets M C O rcg to be dense 
in O reg , if for any x G fl and A G O rcg , x G A there is an N G M such that 
x G iV C A. If, moreover, M is a sublattice such that MM C M holds, we 
call it a dense M -invariant sublattice in O reg . 

Let 9^ C 91 be the image of M through the map Ah- [A]. The following 
fact can be derived as a consequence of density. 

7 Thc latter means i^A]^) = (UA) 1 - =UQUA= HA*. 
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Proposition 8 IfAf is a dense M-invariant sublattice, then the metric space 
(^Mm A r, T Mm A /) is isometric to (Q, d). The isometry is realized by the bijection 
[a] -H- x a . 

Let 91$ C 9\ n C Z(H) be the image of Af through the map Af 3 N h->- 
UN G D\ n . The image is an iM- invariant sublattice in 9t H . The next result 
is a straightforward consequence of the previous one. 

Proposition 9 IfAf C O rcg is a dense M-invariant sublattice then the met- 
ric space (^iM£H#j r iM9t#) ^ s isometric to (Q,d). The isometry is realized by 
the bijection i[a] -H- x a . 

Later, in applications, we will deal with concrete Q and N '. 

The relation (H, d) m = x r iM«^) i s ^ ne fi na l S oa l °f our excursus 

to the lattice theory. It represents the original metric space as collection of 
atoms of relevant Hilbert lattice with inflation. This representation will play 
the key role in reconstruction of Q via inverse data. 

2 Wave spectrum 
2.1 Inflation in % 

10. Basic objects. Let % be a separable Hilbert space, £("H) the lattice 
of its (closed) subspaces equipped with the order ^=C, operations AAB = 
A n B, A V B = {a + b | a G A, b G B}, the complement A A ± = H Q A, 
and extremal elements = {0}, 1 = H. A sublattice of 2(H) is its subset 
closed w.r.t. the operations and complement. Each sublattice contains and 
1. 

By P4 we denote the (orthogonal) projection onto A G £("H). Also, if A 
is a non-closed lineal set, we put P4 := P4. 

Let *B("H) be the algebra of bounded operators. For an S C 23 ("H), by 
ProjiS we denote the set of projections belonging to S. 

For a lattice £ C £("H), with a slight abuse of notation, we put Proj£ := 
{Pa I A G £}. The map £ 3 A > P4 G Proj£ induces the lattice structure 
on Proj£: P A A P B = P AAB , P A V P B = P AvB , (P^ = P M . For P,Q G 
Proj£ the relation P < Q means RanP C RanQ and holds if and only if 
(Px, x) ^ (Qx, x), x EH. The extremal elements of Proj£ are the zero and 
unit operators O and I. 
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The same map relates the order topology on £("%) with the strong oper- 
ator topology on 93 (*H): A = lim A,- in £(%) if and only if P4 = s — IxmP^, 
in <3(H). 

Convention. The metric inflation iM in "H = L 2 ^(fl) is defined on a 
sublattice 9\ n C £("H) (item 8). In contrast to it, in what follows we deal 
with inflations defined on the whole £(H). 

For an inflation / : £(H) -> JFs,{H)i we denote A 1 := (IA)(t), t ^ 0. 
It is also convenient to regard inflation as an operation on projections: / : 

Proj<B(H) ->• Jp r oj<B(«), (IP)(t) = P* := P(/RanP) (i ), t > 0. 

A lattice £ C £(H) is said to be I -invariant if /£ C J-£ holds, i.e. „4 G £ 
implies (IA)(t) G £, t ^ 0. 

Definition 5. Let f C £("H) be a family of subspaces. Define £[J, f] C 
£("H) as the minimal /-invariant lattice, which contains f. 

11. Spectra. Let "H and I be given, £ be an /-invariant lattice. Recall that 
the space of atoms with the 'interaction' topology was introduced in item 3. 

Definition 6. The space fl^. 1 := (fi/£,T/ £ ) is called an atomic spectrum 
of the lattice £. 

There is a version of this notion. Each function / G /£ c Tz is an 
increasing family of subspaces {/(t)}t^o C £("H), , i.e. a nesi pi)]. The 
corresponding nest of projections {Pj} t ^o, Pj '■= Pf(t) determines a self- 
adjoint operator Ef := J^°tdPj. It acts in H and is called an eikonal. The 
set of eikonals is Eik£ := {Ef | / G /£}. 

Definition 7. A metric space f2£ est := {.E^ | a G fij.e} with the distance 
\\E a — Ep\\ is called a nest spectrum of the lattice £. 

Caution! We do not assume E a to be a bounded operators, so that the 
'pathologic' situation dist(E a , Ep) = 00 is not excluded. However, a 'good' 
case, when all the differences E a — Ep are bounded operators, is realized in 
applications. 

One more version is the following. 

Let us say that we deal with a bounded case if the set of eikonals of the 
lattice is uniformly bounded: sup{||i?|| | E G Eik£} < 00. 

With a lattice £ one associates the von Neumann operator algebraic OTc C 
23 ("H) generated by the projections of £, i.e., the minimal von Neumann 
algebra satisfying Proj£ C Proj9t£. 



i.e., a unital weekly closed self-adjoint subalgebra of 03 (Ti): see [T5] , 
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In the bounded case, we have Eik£ C 9t£ (the closure in the strong 
operator topology). The elements of this closure are also called eikonals. 

The set Eik£ is partially ordered: for two eikonals E, E', we write E ^ 
E' if (Ex,x) ^ (E'x,x), x G H. An eikonal E is maximal if E ^ E' implies 
E = E'. By r2g k C Eik£ we denote the set of maximal eikonals. 

Lemma 1 In the bounded case, the set fi| lk is nonempty. 

Proof. By the boundedness, any totally ordered family of eikonals {Ej} has 
an upper bound s-limEj, which is also an eikonal. Hence, the Zorn lemma 
implies fig k ^ 0. □ 

Definition 8. A metric space f2£ lk with the distance \\E — E'\\ is called 
an eikonal spectrum of the lattice £ 

In the general (unbounded) case, one can regularize the eikonals as Ej := 
Jo°° T+it dPf ( e > 0) an d deal with the corresponding spectra fl e ^ k ' e ^ 0. 

Remark. Returning to Definition 6, one more option is to define the 
atomic spectrum as (f2/£, pis) or (Jl/£, aiz) (see sec 1.1, item 3). Our reserve 
of concrete examples is rather poor and provides no preferable choice. 



12. Dynamical system. Let L be a semi-bounded self-adjoint operator in 
H. Without lack of generality, we assume that it is positive definite: 



where h e L x ° c ((0, oo); U) is a "H- valued function of time (control). Its 
solution v = v h (t) is represented by the Duhamel formula 



2.2 Inflation I L 





where dQ\ is the spectral measure of L, x > is a constant. 
Operator L governs the evolution of a dynamical system 

v tt + Lv — h , t > 
v\ t =o = v t \ t =o = 0, 



(2.1) 
(2.2) 



v 




(2.3) 
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(see, e.g., [E])B In system theory, v h is referred to as a trajectory; v h (t) G % 
is a state at the moment t. In applications, v h describes a wave initiated by 
a source h. 

Fix a subspace .A C U. The set := {v h {t) \ h G L l ° c ((0,oo);.A)} of 
all states produced by ^.-valued controls is called reachable (at the moment 
t, from the subspace A). Reachable sets increase: A C B and s ^ t imply 



13. Dynamical inflation. With the system (12.11) . (12.21) one associates a 
map I L : £,(%) -> .F w (/^)(0) := A, (hA)(t) :=V^,t> 0. 

Lemma 2 zs an inflation. 

Proof. The relation (7 L .A)(s) C (I L B)(t) as A C B and < s < t is a 
consequence of the general properties of reachable sets. The only fact we 
need to verify is that the map extends subspaces: A C (I L A)(t), t > 0. 

By X[a,6] we denote the indicator of the segment [a, 6] C R. Fix an r > 
and e G (0,r). Define the functions <^ e (£) := £ -2 X[_ £i£ ](t)sign(— £) and 
y£(t) := <£> £ (t — r + e) for t Gl. Note that f Q r <p r E (t) f (t) dt 4 - /' '(r) as e 
for a smooth /, i.e., <fl(t) converges to S'(t — r) as a distribution. 

For A > 0, define a function 



. f r sm[V\{r-t)] 2 cos(VAe) - cos( v / A2 £ ) - 1 

:= y o ^(tja* = ^ ■ 

Note that if> e (X) —> 1 as £ — > uniformly w.r.t. A in any segment [x, N]. 

e— >0 

Take a nonzero y G A and consider (I2.ip . (12. 2p with the control h e (t) 
(p r £ (t) y. By the properties of ip s one has 



\y - v he (r)\ 



(see (J2H 



, . 00 sin[>/A (r - t)l ,^ r 

y - / / — ^ ^ ^dQ A [<^(%] 

o </x v A 



ipe{X)dQ x y 



[l-V e (A)] dQ A y 



|i-^(A)| 2 d||g A2/ f -> . o. 



9 For x ^ 0, problem (|2.1|) . ()2.2[) is also well defined but the representation (|2.3|) is of 
slightly more complicated form. 
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The order of integration change is easily justified by the Fubini Theorem. 
Thus, y = limv hs (r), whereas v hs (r) G (I L A)(r) holds. Since (I L A)(r) is 

£->0 

closed in H, we get y G (I L A)(r). Hence, A C (I L A)(r), r > 0. □ 

So, each positive definite operator L determines the inflation II, which 
we call a dynamical inflation. 

2.3 Space f^ . 

14. Lattice £^2) and spectra. Fix a subspace 22 G £("H) and say it to be 

a directional subspace. 

Return to the system (I2.ip - fl2.2p . Introduce the class M. v ■ = 
{h G C°° ([0, oo); V) | supp h C (0, oo)} of smooth V- valued controls vanish- 
ing near t = 0. This class determines the sets 



U l v :={ h(t)-v h "(t) 



heMv \ = ( see (E3D ) 



j/i(t) - jf L"5 S in (t- s )L5 



/i G \, t^O, (2.4) 



where ( • )' := 4. These sets are also called reachable. As one can show, the 
sets IA\) increase as t grows. 

Definition 9. The family of subspaces ul,v = {^r)}t>o ^= is called 

a boundary nest. 

The boundary nest determines the lattice £l,v '■= £[Il,Ul,v], which is 
the minimal /^-invariant sublattice in £(H) containing Ul.v (item 10). 

The lattice determines the spectra v and • In the bounded case, 
the spectrum ^ e ^ v is also well defined (item 11). 

15. Lattice £l an d spectra. Let Lq be a closed densely defined sym- 
metric semi-bounded operator with nonzero defect indexes n± = n ^ oo. As 
is easy to see, such an operator is necessarily unbounded. For the sake of sim- 
plicity, we assume it to be positive definite: (L y,y) ^ x|M| 2 , y G DomL 
with x > 0. 

Let L be the Friedrichs extension of Lq, so that L = L* ^ xl and 
L C L C Lq holds j8]. Also, note that 1 ^ dimKerLg = n ^ oo. 

With the operator L one associates two objects: the inflation 1^ and the 
directional subspace T> = KerL^. This pair determines the boundary nest 
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Ul := u L ,KerL ( * = {^kerL h>o and the lattice £ Lo := £ L ,KerL%- The nest and 
lattice determine the corresponding spectra, and we arrive at the key subject 
of the paper. 

Definition 10. The space Qi := fl^ is called a wave spectrum of the 
(symmetric semi-bounded) operator Lq. 

Recall that Ql is endowed with the 'interaction time' topology (item 3). 

By analogy with the latter definition, one can introduce the metric spaces 
Qncst ._ ^nest ^ j n ^ bounded case, Qf k := Q c J k , which are also deter- 

mined by Lq. 

As is evident from their definitions, the spectra are unitary invariants of 
the operator. 

Proposition 10 If U : T-L — >■ H is a unitary operator and L = UL U* then 
Ql o is homeomorphic to f2^ . If H — H 1 © H 2 and L = Lq © Lq then 

^L = ^Lj u ^Lg ■ 

These properties motivate the use of term 'spectrum'. The same properties 
occur for £7^ st and replacing 'homeomorphic' with 'isometric'. 

16. Structures on Ql - The boundary nest Ul can be regarded as an 
element (function) of the space J-z{H)- As such, it can be compared with the 
atoms, which constitute the wave spectrum Ql C F^n)- 

Definition 11. The set 8Ql '■= {a £ I a ^ u io}0 * s sa id to be the 
boundary of Q,l ■ 

Also, it is natural to put dfl L e ^ := dQL - In the bounded case, one 
introduces the boundary eikonal E 9 = J Q °° t dP u t^ ^ and defines d£Vf^ = 

{E G flf-^ | E > E 9 } (see [6]). 

There is a way to represent elements of H as 'functions' on the wave 
spectrum. 

Fix an atom a G VL Lo : a = a(t), t ^ 0. Let P* := P a {t) be the 
corresponding projections. For f,g G "H, we put / = g if there is e = 
e(/, g, a) > such that P*f = P^g as t < e. The relation = is an equivalence. 
The equivalence class [f] a =: G^(a) is called a u>ai>e germ (of the element / 
at the atom a). 
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Recall that a < ul in J-s,(n) means that a(t) C holds for i ^ 0. 
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Definition 12. The germ-valued function G* : a i— >■ [f) a , a G VIl is 
called a wave image of the element /. 

The collection Q := {G^ \ f G %} is a linear space w.r.t. the point-wise 
algebraic operations: (XG? + fiG 9 ) (a) := [Xf + fig] a , a. G f2 io . The linear 
map X : H 3 f > G? E G is called an image operator. 

3 DSBC 

3.1 Green system 

17. Ryzhov's axioms. Consider a collection {H, B; A, T , Vi} of separable 
Hilbert spaces % and B, and densely defined operators A : H — > H and 
r^. : H — )■ (A; = 0, 1) connected via the Green formula 

(Au, v) H - (u, Av) H = (T u, r» B - (Titt, r w) B . 

The space % is called an inner space] B and are referred to as a boundary 
values space and the boundary operators respectively [H] . Such a collection 
is said to be a Green system. 

The following additional conditions are imposed. 
Rl. DomTfc D Domi holds. The restriction AlxerronKeiT! =: L is a densely 
defined symmetric positive definite operator with nonzero defect indexes. 
The relation A = L* Q is valid ('bar' is the operator closure). 

R2. The restriction A|Kerr =: L coincides with the Friedrichs extension of 
L , so that we have L C L C Lq = A. Operator L~ l is bounded and defined 
on H. 

R3. The subspaces A := Ker A and V := Ker Lq are such that the relations 
A = V and IV4 = B hold. 

These conditions were introduced by V.A.Ryzhov [IB] , which puts them 
as basic axioms. Note, that there are a few versions of such an axiomatics 
but the one proposed in [12] is most relevant for applications to forward and 
inverse multidimensional problems of mathematical physics. 

The following consequences are derived from Rl-4 [IB] . 
CI. The operator II := (TiL^ 1 )* : B — > H is bounded. The set Ran II is 
dense in T>. 

C2. The representation A = {y G Domi | HT y = y} is valid. 
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C3. Since L is the extension of Lq by Friedrichs, the relations DomLo = 
L~ l \H V] and L = L\L-i[ueT>] easily follows from the definition of such 
an extension (see [8]). 

18. Illustration. Let Q be a C°°-smooth compact Riemannian man- 
ifold with the boundary T, A the (scalar) Beltrami-Laplace operator in 
H := L 2 (fi), v the outward normal on T, B := L 2 (T). Denote l"1 A = 
— A\ H 2rQ\, T u = u\r, Tiu = d u u\r , so that r 0il are the trace operators. 
The collection {H, B; A, T , Ti} is a Green system. Other operators, which 
enter in Ryzhov's axiomatics, are the following: 

• L = — A\ H 2^ is the minimal Laplacian that coincides with the closure of 
—A\c^(n\dn) 

• L = — A| 

H 2 (n)nH^(n) 1S the self-adjoint Dirichlet Laplacian 

• Lq — — A|{ yG ^ | Ayeu} is the maximal Laplacian 

• A = {y G H 2 (Q) | Ay = 0} is the set of harmonic functions of the class 

• T> = {y G % | Ay = 0} is the subspace of all harmonic functions in L 2 {Vt) 

• II : B — > H is the harmonic continuation operator (the Dirichlet problem 
solver): Hip = u is equivalent to Au = in Q, u\r = (p. 

3.2 Evolutionary DSBC 

19. Dynamical system. The Green system determines an evolutionary 
dynamical system with boundary control 

u a + Au = in U, < t < oo (3.1) 

u\ t =o = u t \ t =o = in H (3.2) 

T u = f{t) in B, «C t < oo, (3.3) 

where / is a boundary control, u = u^(t) is the solution (wave). The space 
of controls J 7 = L^ ((0, oo); B) is said to be outer. 

Assign / to a class J- + C T if it belongs to C°° ([0, oo); B), takes the values 
in r Domv4 C B, and vanishes near t = 0, i.e., satisfies supp / C (0,oo). 
Also, note that / G J-+ implies II (/( • )) G M-d (see item 14). 

ll H k are the Sobolev classes; H%{Q) = {y e H 2 {VL)\y = \Vy\ = onT}; d v is the 
differentiation w.r.t. the outward normal on T. 
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Lemma 3 For f G J-+, the classical solution u* to problem (13. lft — H 3 . 3 1) is 

represented in the form 



u f(t) = h{t) 



L 2 sin 



{t-s)Li 



h"(s) ds. 



t > 



(3.4) 



with h := n (/(•)) G -Mo. 



Proof. Introducing a new unknown w = w^(t) := u*{t)—Tl [f{t)) and taking 
into account CI (item 17), we easily get the system 



w tt + Aw = -U{f tt (t)) 
w\ t =o = w t \t=o = 

r w = o 



in Ji, < t < oo 
in U 

in B, < t < oo . 



With regard to the definition of the operator L (see the axiom R2), this 
problem can be rewritten in the form 



w tt + Lw = -h tt 
w\ t=0 = w t \ t=0 = 

and then solved by the Duhamel formula 



in H 1 < t < oo 
in U 



w f (t) 



L 2 s i 



sin 



(t-s)Li 



h"(s) ds . 



Returning back to w = w' + II/, we arrive at (13.41) . □ 
20. Reachable sets. The sets 
U\ := {u f {t) | / G J^+j = ( see (E3D ) = 



h(t) 



L 2 sin 



(t-s)L3 



h"(s) ds 



h = Uf(-), fe^ 



t ^ o 

(3.5) 



are said to be reachable from boundary. 

The Green system, which governs the DSBC, determines the certain pair 
L,V, which in turn determines the family {U^} by (12. 4p . Comparing (12.4D 
with (13. 5p . we easily conclude that the embedding IA 1 + C Uj, holds. Moreover, 
the density properties R3 (item 17) enable one to derive U\_ = U^,, t ^ 0. 
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It is the latter relation, which inspires the definition (12 Ah and motivates the 
terms 'reachable sets', 'boundary nest', etc in the general case (item 14), 
where neither boundary value space nor boundary operators are defined. 

21. Illustration. Return to the item 19. The DSBC (j3TT]) -(l3~15 ]) associ- 
ated with the Riemannian manifold is governed by the wave equation and is 
of the form 

u tt - Am = in x (0, oo) (3.6) 

u\ t=0 = u t \ t =o = in Q (3.7) 

u|r = /(*) for ^ t < oo (3.8) 

with a boundary control / G T = L 1 ^ ((0, oo); L 2 (T)). The solution u = 
uf(x,t) describes a wave, which is initiated by boundary sources and propa- 
gates from the boundary into the manifold with the speed 1. For / G J-+ = 
C°° ([0, oo); C°°(r)) provided supp/ G (0, oo), the solution is classical. 

By the finiteness of the wave propagation speed, at a moment t the 
waves fill a near-boundary subdomain r* := {x G Q | dist (x,T) < t}. Cor- 
respondingly, the reachable sets U\_ increase as t grows and the relation 
U\ C UY\ t ^ holds Closing in H, we get C HT\ t > 

So, if the pair L, V (or, equivalently, the operator L ) appears in the 
framework of a Green system, then introduced by the general definition 
(12. 4p can be imagine as the sets of waves produced by boundary controls. The 
question arises: What is the meaning of the corresponding wave spectrum 
^l,x> (= ^L )? m a sense, it is the question, which this paper is written for. 
The answer (section 3) is that, in generic cases, Ql is identical to Q. 

22. Boundary controllability. Return to the abstract DSBC (13. ip ~ 
(13. 3 p and define its certain property. Begin with the following observation. 
Since the class of controls T+ satisfies -^■J r + = the reachable sets (13. 5p 
satisfy AU\_ = U\. Indeed, taking / G J-+ we have 

Au f (t) = ( see $£T$ ) = -u{ t {t) = u~ f " \t) G U\ . (3.9) 

By the same relations, u/(t) = Au 9 (t) holds with g = — (J *) 2 / G Hence, 
the sets reduce the operator A and its parts v4| w t are well defined. 



12 Geometric subspaces HA are defined in item 8. 
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Definition 13. The DSBC (I3.1l) -( l3~3l) is said to be controllable from 
boundary at the time t — T if A\ u t = A holds, i.e., one haJ^I 

{{uf(T),Auf(T)} | / EJ r +} = graphs = (see Rl) = graph L*. (3.10) 

Controllability means two things. First, since A is densely defined in "H, 
the equality (I3.10P implies U\_ = H, t ^ T , i.e., for large times the reachable 
sets become rich enough (dense in H). Second, the 'wave part' A\ u t of 
the operator A, which governs the evolution of the system, represents the 
operator in substantial. 

In applications to problems in bounded domains, such a property 'ever 
holds' (typically, for large enough times T). In particular, the system (13. 6ft - 
(13.81) is controllable from boundary for any T > maxdist (x, T) [2], [I]. 

Let us represent the property (I3.10p in the form appropriate for what 
follows. 

Restrict the system (I3.ip - (I3.3I) on a finite time interval [0, T\. Define the 
Hilbert space of controls J-" T = L 2 ([0,T];£>) and the corresponding smooth 
class J 7 ^ C J-" T . 

Introduce a control operator W T : J-" T — > H, DomW T = J-T, W T f : = 

v,f(T). Let W T = U T \ W T \ be its polar decomposition, where \W T \ : = 
i 

((W T )*W T ) 2 acts in J rT , and U T is an isometry from Ran |iy T | C T T onto 
RanI^ T C U (see, e.g., [8]). 

Lemma 4 // the DSBC (13. ip — (13. 3 p is controllable att = T then the relation 
{{\W?\f, \WT\{-f»)} | feF + } = (U t )*LqU t holds. 

Proof. Represent (1 3J0J) in the equivalent form {{W T f, W T (-f")} \ f G 7+} 
= graph Lq. Since Ran?7 T = U T = H, the isometry U T is a unitary oper- 
ator. Applying it to the latter representation, one gets the assertion of the 
lemma. □ 

As a consequence, we conclude the following. 

Proposition 11 If the DSBC (13. ip — (13. 3p is controllable at t = T then the 
operator \ W T \ determines the operator Lq up to unitary equivalence. 

23. Response operator. In the DSBC ([3l])-([33D restricted on [0,T], an 
'input-output' correspondence is described by the response operator R T : 
J= T Domi? = 1%, (R T f) (t) := T 1 (u f (t)) , < t ^ T. 

13 Below the closure is taken va.T-Ly.T-L] graphs := {{y, Ay} \ y £ DomA}. 
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Illustration. The response operator of the DSBC f l3.6p - fl3.8l) is R T : / t— > 

dvw\rx[Q,T\- 

The key fact of the BC-method is that the operator R 2T determines the 
operator C T := (W T )*W T through an explicit formula [2],[3],[I]. 

Proposition 12 The representation C T = ^(S T )* R 2T J 2T S T holds, where 
the operator S T : T T — > T 2T extends controls from [0, T] to [0, 2T] by oddness 
w.r.t. t = T, J 2T : 7 2T ->• ^ 2T , (J 2T /)(t) = /„' /(s) ds. 

Hence, R 2T determines the modulus \W T \ = (C T )^ . By Proposition [TT| we 
conclude that R 2T determines the operator Lq up to unitary equivalence. 
Since Lq = Lq*, we arrive at the following basic fact. 

Proposition 13 If the DSBC ( 13. lft — H 3 . 3 1) is controllable from boundary at 
t = T then its response operator R 2T determines the operator Lq up to unitary 
equivalence. 

24. Illustration. The system ( I3.6p - fl3.8p is also controllable from boundary. 
Such a property is a partial case of the following general fact. 

Return to the system ( I2.ip -( 12TT]1 . In our case, the operator L governing 
its evolution is the Dirichlet Laplacian —A (item 18). Fix a set A e O reg . 
The reachable sets V l nA consist of the waves produced by sources supported 
in A C Q. Since the waves propagate with unit velocity, the embedding 
Vha — HA 1 holds evidently. The character of this embedding is a subject of 
control theory of hyperbolic PDE. 

The principal result is that the relation V l UA = HA 1 is valid for any 
A G O veg and t ^ 0. It is derived from the fundamental Holmgren- John- 
Tataru uniqueness theorem (see, e.g., [2], jl]). In control theory this property 
is referred to as a local controllability of manifolds. In notation of item 13, 
it takes the form: (I L HA)(t) = HA* holds for any A e C reg , t ^ 0. Since 
T-LA 1 = (iM7iA)(t) by the definition of metric inflation on £H W (item 8), we 
arrive at the following formulation of the local controllability. 

Proposition 14 The inflations Ii and iM coincide on the lattice 

Return to the system fl3.6p - fl3.8p and the embedding U^, C T-LY 1 (item 
21). The same HJT-theorem implies the equality ICp = HT*, t ^ 0, which is 
referred to as a local boundary controllability of the manifold Q. 
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Recall that the boundary nest ul = {Wp}t>o (T> = KerLg) is defined 
in item 15. Let b = {r*}^ C O reg be the family of metric neighborhoods 
of the boundary Y. Denote [b] = {[r*]}^o C £H (items 7,8). Boundary 
controllability of Q is equivalent to the following. 

Proposition 15 The relation i[T T ] = U^, t ^ holds. Hence, i[b] = Uc . 

Boundary controllability implies the following. Since the family {T*} 
exhausts fl for any T ^ T* := sup,j. gf2 d(x, T), the boundary nest {U^} t ^T 
exhausts the space Tl as T ^ T*. By this, the system (I3.6p - (l3.8p turns out 
to be controllable as T ^ T* [2], [I]. 

Hence, by Proposition [T3| given for a fixed T 2T* the response operator 
R T of the system fl3.6p - fl3.8p determines the minimal Laplacian L up to 
unitary equivalence. 



25. Weyl function. Here we follow the paper [16], and deal with the same 
Green system {H, B; A, T , Ti} and the associated operators L , L (item 17). 

The problem 



is referred to as a stationary DSBC. For ip 6 ToDomA and z G C\specL, 
such a problem has a unique solution u = u^, which is a Dom^ -valued 
function of z. 

The 'input-output' correspondence in the system fl3.lip - fl3.12p is realized 
by an operator-valued function W(z) : B — >■ B, W(z){p := Tiuf (z ^ specL). 
It is called the Weyl function and plays the role of data in frequency domain 
inverse problems. 

The following important fact is established in [16]. Recall that a sym- 
metric operator in T-L is said to be completely non- self adjoint if there is no 
subspace in "H, in which the operator induces a self-adjoint part. 

Proposition 16 If the Green system is such that the operator Lq is com- 
pletely non- self adjoint, then the Weyl function determines the operator L 
up to unitary equivalence. 



3.3 Stationary DSBC 



(A - zl) u = 

r u = ip 



in ft, zeC 
in B 



(3.11) 
(3.12) 



23 



26. Illustration. Return to item 17. The DSBC f l3TTTD -f l3TT2l) associated 
with the Riemannian manifold is 

(A + z)u = mQ (3.13) 
u|r = ¥>, (3-14) 

where A = -A|#2 (c) . 

Lemma 5 The operator L = — A\ H 2^ is completely non-selfadjoint. 

Proof. Assume that there exists a subspace K, C H such that the operator 
Lq := — A|^ nif 2^ 7^ O is self-adjoint in /C. In the mean time, Lq is a part 
of L, which is a self-adjoint operator with the discrete spectrum. Hence, 
specL^ is also discrete; each of its eigenfunctions satisfies — A0 = A0 in Q 
and belongs to Hq(Q). The latter implies = d u (p = on T. This leads to 
<p = by the well-known E.Landis uniqueness theorem for solutions to the 
Cauchy problem for elliptic equations. Hence, Lq = O in contradiction to 
the assumption. □ 

The Weyl function of the system is W(z){p = d u u^\r (z ^ specL). By 
the aforesaid, the function W determines the minimal Laplacian Lq of the 
manifold Q up to unitary equivalence. 

27. Spectral data. Besides the Weyl function, there is one more kind 
of boundary inverse boundary data associated with the DSBC (I3.13p -f l3.14p . 
Let {Afc}^ : < Ai < A2 ^ A3 ^ • ■ ■ — > 00 be the spectrum of the Dirichlet 
Laplacian L. Let {(p^kLi '■ L(pk = ^k4>k be the corresponding eigen basis in 
U normalized by (<p k , <pi) = 5 kl . 

The set of pairs := {A^; ct^lr}^ is called the (Dirichlet) spectral 
data of the manifold Q. 

The well-known fact is that these data determine the Weyl function and 
vice versa (see, e.g., [IS]). Hence, determines the minimal Laplacian L 
up to unitary equivalence. However, such a determination can be realized 
not through W but in more explicit way. 

Namely, let U : H — > H := h, Uy = y := {{y , 4>k)}kLi be the Fourier 
transform that diagonalizes L: L := ULU* = diag{Ai, A2 , . . . }. For any 
harmonic function a G A, its coefficients are (a, 4>k) = — j- J T a dv4>k dT 
that can be verified by integration by parts. Therefore, the spectral data En 

determine the image A := UA C H and its closure V = UV = A. Thus, the 
determination =^> L , T> occurs. 
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In the mean time, the relation C3 (item 17) implies Lq = U*LqII = 
L\z- 1 [neT>] isometry of U. Thus, L is a unitary copy of L constructed 
via the spectral data. 

4 Reconstruction of manifolds 

4.1 Inverse problems 

28. Setup. In inverse problems (IP) for DSBC associated with manifolds, 
one needs to recover the manifold via its boundary inverse datep Namely, 
IP 1. given for a fixed T > 2maxdist (x, T) the response operator R T of the 
system fl3.6l) - fl3.8p . to recover the manifold Q 

IP 2. given the Weyl function W of the system (I3.13I) - (I3.14I) . to recover the 
manifold Q 

IP 3. given the spectral data En, to recover the manifold fl 

The problems are called time- domain, frequency- domain, and spectral respec- 
tively 

Setting the goal to determine an unknown manifold from its boundary 
inverse data, we have to keep in mind the evident nonuiqueness of such a 
determination: all isometric manifolds with the mutual boundary have the 
same data. Therefore, the only reasonable understanding of 'to recover' is to 
construct a manifold, which possesses the prescribed data [I]. 

As we saw, the common feature of problems IP 1-3 is that their data 
determine the minimal Laplacian L up to unitary equivalence. By this, 
each kind of data determines the wave spectrum Ql up to isometry (see 
Proposition [TO]) . As will be shown, for a wide class of manifolds the relation 

l = ^ holds. Hence, for such manifolds, to solve the IPs it suffices to 

extract a unitary copy L from the data, find its wave spectrum f2^ o fi io , 
and thus get an isometric copy of Q. It is the program for the rest of the 
paper. 

29. Simple manifolds. Recall that we deal with a compact smooth 
Riemannian manifold Q with the boundary T. The family b = {r*}t>o con- 
sists of metric neighborhoods of T. Nets and dense lattices were introduced 

14 In concrete applications (acoustics, geophysics, electrodynamics, etc), these data for- 
malize the measurements implemented at the boundary. 



25 



in item 9. £[M, b] C O ieg is the minimal M-invariant (sub)lattice, which 
contains b. 

We say fl to be a simple manifold if the lattice £[M, b] is dense in O reg . 

The evident obstacle for a manifold to be simple is its symmetries. For a 
ball Q = {x G M ra | |x| ^ 1}, the lattice £[b, M] consists of sums of 'annuluses' 
of the form {x G Q | ^ a < \x\ < b ^ 1}. Surely, such a system is not a 
net in the ball. A plane triangle is simple if and only if its legs are pair-wise 
nonequal. Easily checkable sufficient conditions on the shape of Q C M. n , 
which provide the simplicity, are proposed in pQ. They are also appropriate 
for Riemannian manifolds and show that simplicity is a generic property: it 
can be provided by arbitrarily small smooth variations of the boundary T 0. 



30. Solving IPs. The following result provides reconstruction of f2. 

Theorem 1 Let Q be a simple manifold, L = — A\ H 2^ the minimal Lapla- 
cian, Ql its wave spectrum. There exists an isometry ( of metric spaces ) i* 
that maps Ql onto Q, the relation z*(<9fi£ ) = T being valid. 

Proof. Denote [b] := {[r*]} t>0 C 91 Let £[M, [b]] C 91 be the image of 
£[M, b] through the 'projection' A i-> [A] (item 7). 

Propositions MEB imply i£ [M, [b]] = £[zM,z[b]] = £[/l,u Lo ] = £ Lo C 

Taking into account the simplicity condition and applying Proposition [9] 
to the case Af = £[M, b], we conclude that Ql is isometric to (f2,d). The 
isometry is realized by the bijection i* : i[a] x a . 

To compare the atoms i[a], which constitute f2i , with the boundary nest 
Ul is in fact to compare the metric neighborhoods {x a Y with the metric 
neighborhoods {T*}. Since {x a Y C P, t ^ is valid if and only if x a G T, 
we conclude that i*(<9f2 Lo ) = V. □ 

Thus, to solve the IPs 1-3 in the case of simple Q, it suffices to determine 
(from the inverse data) a relevant unitary copy L of the minimal Laplacian, 
and then find its wave spectrum fi^ o . 

31. Remarks. 

15 Presumably, any compact manifold with trivial symmetry group is simple but it is a 
conjecture. In the mean time, for noncompact manifolds this is not true. 
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1. Regarding non-simple manifolds, note the following. If the symmetry 
group of Q is nontrivial then, presumably, Ql q is isometric to the properly 
metrized set of the group orbits. Such a conjecture is motivated by the 
following easily verifiable examples. 

• For a ball Q = {x G M™ | \x\ ^ r}, the spectrum Ql is isometric to the 
segment [0,r] C R. Its boundary 8Ql is identical to the endpoint {0}. 

2 2 

• For an ellipse = {(x, y) G IR 2 | ^2 +f? ^ 1}, * s isometric to its quarter 

2 2 

Q PI {(x,y) I x ^ 0, y ^ 0}, whereas <9f2L is isometric to {(x,?/) | ^ + = 
1, x ^ 0, y ^ 0}. 

• Let to C {(xi,X2) G M 2 I X2 > 0} be a compact domain with the smooth 
boundary. Let Q be a torus in M 3 , which appears as result of rotation of to 
around the xi-axis. Then VLl to and 0^1l =" 9a;. 

2. In applications a possible lack of simplicity is not an obstacle for solv- 
ing problems IP 1-3 because their data not only determine a copy of L 
but contain substantially more information about Q. Roughly speaking, the 
matter is as follows. When we deal with these problems, the boundary T is 
given. By this, instead of the boundary nest Ul of the sets reachable from 
the whole T (see (13. 5J) ). we can use the much richer family u' Lq = {Wj}^o,ucr 
of sets reachable from any patch a C T of positive measure 0. Therefore, 
even though the density of the lattice £[Il, Ul ] in may be violated by 
symmetries, the lattice £[Il, u' l J is always dense. As a result, the wave spec- 
trum corresponding to the dense lattice turns out to be isometric to Q. The 
latter is the key fact, which enables one to reconstruct Q: see [5] for detail. 

3. The spectra f^ st and flf^ are also appropriate for reconstruction. If 
n is simple, one has Vl Lo is = fi£f is = m fif* is = m (fi, d) [5], [6]. 

4. If Q is noncompact, the definition of simplicity remains to be mean- 
ingful, local controllability is in force, and H = U t> oU^, holds. One can show 
that the response operator R T known for all T > determines the simple 
manifold up to isometry. Also, defining mutatis mutandis the Weyl function 
and spectral data for a noncompact Q, one can obtain the same result: these 
data determine the simple manifold up to isometry. 

32. Algebras in reconstruction. Recall that the von Neumann algebra 
91c C 53(H) associated with the lattice £ C £{%) was introduced in item 

16 More precisely, consists of the solutions (waves) u^(t) produced by the boundary 
controls / supported on a x [0, 00) 
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11. In the bounded case, along with OTo one can define the algebra as 
the minimal norm-closed subalgebra of 03(7/), which contains all maximal 
eikonals. 

For the algebras = 91e Lo and <Tl ; = £z Lq associated with manifold, 
the following holds [5], [6]. 

(i) Both of these algebras are commutative. The embedding Cl C OT^o is 
dense in the strong operator topology in 03 (%). 

(ii) If Q is simple then €l is isometrically isomorphic to the algebra C(Q) 
of continuous functions. By this, its spectrurrf^l <£x is homeomorphic to fl 

These results are applied to reconstruction by the scheme {inverse data} =>■ 

Note that commutativity is derived from local controllability of the sys- 
tem (I3.6p - fl3.8p . In the corresponding dynamical system on a graph, a lack of 
controllability occurs and, as a result, these algebras turn out to be noncom- 
mutative 0. This leads to problems and difficulties in reconstruction, which 
are not overcome yet. In particular, the relations between the spectra Qi 
and <£l are not clear. 

4.2 Comments 

33. A look at isospectrality Let spec L = {Xk}^ =1 be the spectrum of the 
Dirichlet Laplacian on Q (item 27). The question: "Does specL determine 
fl up to isometry?" is a version of the classical M.Kac's drum problem[12|. 
The negative answer is well known (see, e.g., [9]) but, as far as we know, 
the satisfactory description of the set of isospectral manifolds is not obtained 
yet. The following is some observations in concern with such a description. 

Assume that we deal with a simple fl In accordance with Theorem 1, such 
a manifold is determined by any unitary copy Lq of the operator Lq C L. If 
the spectrum of L is given, to get such a copy it suffices to possess the Fourier 
image T> = UT> of the harmonic subspace in Ti = h- see C3, item 170 In the 
mean time, as is evident, if Q and Q' are isometric, then the corresponding 
images are identical: T> = V . Therefore, Q and Q' are isospectral but not 
isometric if and only if T> ^ T>' . In other words, the subspace T> is a relevant 
'index', which distinguishes the isospectral manifolds. 

17 i.e., the set of maximal ideals of €l [15] . 
18 N.Wada, private communication. 
19 It is the fact, which is exploited in 1 
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As an image of harmonic functions, which is admissible for the given L ■ 
diag {Ai, A2, • • • }, a subspace V C h has to obey the following conditions: 



A lineal set £^ := L 1 



is dense in l 2 , whereas replacement of D 

by any wider subspace V D V leads to the lack of density: clos 7^ h, 

2. Extending an operator L\c by Friedrichs, one gets L. 

In the mean time, taking any subspace V <Z l 2 obeying 1,2 0, one can 
construct a symmetric operator L by C3, and then find its wave spectrum 
as a candidate to be a drum. However, the open question is whether 
such a 'drum' is human (is a manifold). 

34. Wave model. Return to the abstract system ( I3.ip -( 13~U]) and assume it 
to be controllable at t = T. Reduce the system to the interval ^ t ^ T. 
Recall that the image and control operators X : "H — > Q and W T : T T — > Ti 
were introduced in items 16 and 22 respectively. The composition V T := 
XW T : J-" T — > Q is called a visualizing operator [2], [3], [1]. 

Let the response operator R 2T be given. The following is a way to con- 
struct a canonical 'functional' model of the operator Lq. 

• R 2T determines the operator \W T \ in J-" T (item 23). In what follows, it is 
regarded as a model control operator W T := |M^ T |, which acts from to a 
model inner space H := T T . 

• Determine the operator Lq in Iri as the operator of the graph 
|{H> T /, W T (-f")} I / G 7^+} (Lemma H item 22). Find L = Lq*. 

• Find the wave spectrum Qi Q and recover the germ space Q on it. Determine 
the image operator X : 7-L — > Q. Compose the visualizing operator V T = 

iw T -.T T ->g. 

• Define (L nod )* as an operator in Q determined by the graph 

{{V T f, V T {-f")}\feF + ). 

Surely, it is just a draf10of the model and plan for future work: one needs 
to endow the germ space Q with relevant Hilbert space attributes. Presum- 
ably, in 'good cases', Q = L 2iAt (fiL ). Also, the model operator is expected 
to be local: supp(L nod )*?/ C suppy, whereas the model trace operators r 0j i 



20 such subspaces do exist (M.M.Faddeev, private communication) 
21 Some detail see in [6], sec 3.4 . 
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are connected with the restriction y i-)- y|an Lfl S Hopefully, the collection 
{Q, B; (L™ od )*, r , r\} constitutes the Green system, which is a canonical 
model of the original {H, £>; A, F Q , Ti}. The model is determined by R 2T . 

Such a model is in the spirit of general system theory [13] , where it would 
be regarded as a realization relevant to the transfer operator function R 2T . 
Remarkable point is the role of a time in its construction. 

35. Open question. For any operator L of the class under consideration, 
the lattice £l is a well-defined object, £l ^ {0} being hold. We have 
neither a proof nor a counterexample to the following principal conjecture: 

7^ 0. However, there is example of the operator L such that VLl consists 
of a single point. 

36. A bit of philosophy. In applications, the external observer pursues 
the goal to recover a manifold Q via measurements at its boundary T. The 
observer prospects Q with waves u* produced by boundary controls. These 
waves propagate into the manifold, interact with its inner structure and 
accumulate information about the latter. The result of interaction is also 
recorded at T. The observer has to extract the information about Q from 
the recorded. 

By the rule of game in IPs, the manifold itself is invisible (unreachable) in 
principle. Therefore, the only thing the observer can hope for, is to construct 
from the measurements an image of Q possibly resembling the original. By 
the same rule, the only admissible material for constructing is the waves w . 
To be properly formalized, such a look at the problem needs two things: 

• an object that codes exhausting information about Q and, in the mean 
time, is determined by the measurements 

• a mechanism that decodes this information. 

Resuming our paper, the first is the minimal Laplacian Lq, whereas to decode 
information is to determine its wave spectrum constructed from the waves 
. It is n^g, which is a relevant image of Q. 

The given paper promotes an algebraic trend in the BC-method [5], by 
which to solve IPs is to find spectra of relevant lattices and algebras. An 
attempt to apply this philosophy to solving new problems would be quite 
reasonable. An encouraging fact is that in all above-mentioned unsolved 



22 As far as we know, the known models of symmetric operators do not possess such 
properties [17] ■ 
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IPs of anisotropic elasticity and electrodynamics, graphs, etc, the wave spec- 
trum Ql does exist. However, to recognize how it looks like and verify (if 
true!) that Ql is isometric (homeomorphic) to Q is difficult in view of very 
complicated structure of the corresponding reachable sets U l . 

4.3 Appendix 

37. Basic lemma. Recall the notation: for a set A C f2, A is its metric 
closure, mtA is the set of interior points, A 1 is the metric neighborhood of 
radius t, A := A. If A e O then A C mtA and A = mtA holds. 

Return to item 5. Let / = f(t),t ^ be an element of MO. Define 
the set / := n t>0 f(t) C fi. Define the functions /*(£) = (Mf)(t) = /* as 

t > 0, /*(0) = /(0) and f*(t) = mtp, t > 0. 

Lemma 6 (i) If f ^ 0j- o t/ien / = / ^ and tae relations /* ^ / ^ /* 
/ioW m J"e>. (uj /// and g satisfy f = g then f(t) = g(t) (= /*) as t ^ 0. 

Proof. 

1. If / \ /j- G MO then f(t) = U m fj(t), t ^ 0. Therefore, f k (0) C 
U^i/,(0) C U i>1 (/ i (0)) t = U 3>l f j (t) = f(t), t > 0. Hence, / D / fc (0) ^ 0. 

2. If / / fj e MO then f(t) = int n 3 >i fj(t), t ^ 0. Define a closed set 
F = n m fj(0) C n and show that F ^ 0. 

Assume F = 0. Since / J+ i(0) C /j(0), for any x G Q and t > there 
is Jo = jo(x,t) such that {a;}' fl /j(0) as j > jo- Indeed, otherwise, by 
assumptions Al,2, the ball {x} 1 has to contain the points of F. Hence, 
x ^ (/j(0))* as j > jo- Since x is arbitrary, we have = Uj^i(/j(0))* = 
Uj>i(fj{0)Y = U m fj{t). Therefore /(*) = int n m fj(t) = $,\e., f{t) = 
0o, t ^ 0. It means that / = 0jr o in contradiction with assumptions of the 
lemma. So, F ^ 0. 

3. Show that F — f, i.e., F does not depend on {/j)n For every j ^ 1, 
we have J^O) = n t >o(/j(0))* = fW,(*) 2 fW(*) = /■ Hence F = 
' ; !/,((') 2f- 

On the other hand, the monotonicity fj + i(0) C fj(0) implies 
F = n m fj(0) C (n^i/j^O))* C H^i (/,-(()))'. Since the next to the last set 
is open as t > 0, we have F C (n J - > i/ J -(0))* C int D 3 >i (fj(0)Y Q int H^i 



However, the limit / can depend on {fj}: there are examples for fl = M"! 
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(ft(0)Y C int n m fj(t) = fit) for all t > 0. Hence F C n t>0 /j(t) = /, and 
we arrive at F — f. 

Thus, we obtain F = f ^ 0. 

4. Show that /* < /. Choosing MO 3 ft \ /, for t > one has /* = 
C (f0jy = (f 3 (0)Y = ft(t). This implies /* C D m ft(t). Since f is 

an open set, the embedding /* C int U,>i /,(£) = /(£) holds. The latter 
means that ft(t) < /(t) in (9 as t > 0. The definition of /* at t — leads to 
/„(t) < /(t), i > in O, i.e., /* ^ / in T Q . 

Show that / ^ /*. Choose MO 3 ft \ f that means /(*) = int n m 
ft(t), t ^ 0. For t = one has /(0) = int n^i £(0) C int n^i/^O) = int/ = 
int/ = /*(0). For t > 0, with regard to monotonicity of {/j}4, we have 

/(t)_= intrw^t) = intn^x (fMY = int (^(O))* c int(n^5fo^ = 

int/* = /*(*). Hence / ^ /* is valid. 

Thus, the part (i) of the lemma is proven. 

5. For t > 0, since /* is an open set, one has /* = int/*. Therefore, f*(t) = 

/*(*) = /*, and (i) implies 7(7) = /*. Hence, J(t) = f = J 1 = gjt) as * > 0. 
Let t = 0. Choosing MO 3 ft \ /, one has /(0) = int n^i ft(0) C 

int /j(0) = int/. Hence 7(0) C int/ C /. Show that 7(0) = int/. 
Indeed, assuming the opposite, one can find x G / separated from /(0) with a 
positive distance. In the mean time, defining f £ by f £ (t) = (/(0)) e+t , t ^ 0, 
we get /*(£) = /* C /(t) c / e (t). However, the relation /* C / e (t) is 

impossible for small enough t and e by the choice of x. Hence, /(0) = int/ 

does hold. 

The latter implies /(0) = int/ = mtg = g(0). Thus, we get f(t) = g(t) 
for all t ^ and prove (ii). □ 

Key words: symmetric semi-bounded operator, lattice with inflation, evo- 
lutionary dynamical system, wave spectrum, reconstruction of manifolds 
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